arXiv:1508.02673vl [math.AP] 11 Aug 2015 
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Abstract. Optimal partial transport, which was initially studied by 
Caffarelli and McCann [3], is a variant of optimal transport theory, where 
only a portion of mass is to be transported in an efficient way. Free 
boundaries naturally arise as the boundary of the region where the actual 
transport occurs. This paper considers the evolution dynamics of the 
free boundaries in terms of the change of m, the allowed amount of 
transported mass or the change of A, the transportation cost cap, i.e. the 
allowed maximum cost for a unit mass to be transported. Focusing on 
the quadratic cost function, we show Holder and Lipschitz estimates on 
the speed of the free boundary motion in terms of m and A, respectively. 
It is also shown that the parameter m is a Lipschitz function of A, which 
previously was known only to be a continuous increasing function [3]. 


1. Introduction 

Given two mass distributions we consider the phenomena of matching 
them together in a cost efficient way. We are interested in the case when 
only a fraction of the mass is to be matched and therefore only a part of the 
mass distributions are transported. This, so-called, optimal partial transport 
problem, has been an interest among researchers starting from the work of 
Caffarelli and McCann [3], It is a natural generalization of the optimal 
transport problem of Monge and Kantorovich mm where the full masses 
are matched; see, P21EL9] for a modern survey. 

In this paper, our aim is to understand the dynamical behaviour of the 
solution to the optimal partial transport problem: First, one can formulate 
this partial transport problem in two equivalent forms, using two different 
parameters. One is the amount m of mass to be transported, and the other 
one is the transportation cost cap A, the allowed maximum cost for a unit 
mass to be transported. See Section [2] for more details for the two equivalent 
formulations of the partial transport problem. As these parameters increase, 
the active region where the transport of mass actually occurs, changes (in 
fact, it changes monotonically 0 ). Focusing on the cost function c(x , y ) = 
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\\x — y | 2 , the main goal of this paper is to present certain quantitative 
estimates (see Theorem 14.51 for precise statements) on the change of the 
active region with respect to the change of the parameters, m or A. In fact, 
such an estimate with respect to m is not difficult, since the difference in 
amount of mass can be easily related to the volume change of the active 
region. Estimates with respect to A are more difficult to obtain, since the 
relation between the transportation cost cap and the amount of transported 
mass, thus with the volume of active region, is indirect. Note that the 
parameter m can be considered as a function of A, and it is known to be 
monotone 0. As a byproduct of our estimates, we show that A i->- m( A) is 
Lipschitz (see Theorem 14.51 item 4, and see also, Corollary 14.61) . Our main 
technical tool is a certain monotonicity of the potential functions associated 
to the partial transport problem (see Theorem l6.il) . 

Even though we do not pursue it here, these results can be extended 
to more general class of examples (see Remark 13.2|) . However, to show a 
strict separation of the free boundaries for different values of parameters, 
we use the special structure of c(x,y) = \\x — y | 2 , and the assumption 
that mass distributions are smooth: we show a version of strong maximum 
principle for the Monge-Ampere equation, from which we show that the ac¬ 
tive region grows strictly monotonically in a point-wise sense (see Section [8j 
Theorem 18.5p . 


Organization of the paper: Section [2] gives the notation throughout the 
paper, and explains the equivalence between the two approaches to partial 
transport problems, with respect to m and A, respectively. Section [3] sets up 
the main assumptions, and explains preliminary results. Section [4] explains 
the main results, whose proofs are given in Sections [5] and [71 where the 
key monotonicity of potential functions is given in Section [6) Finally, Sec¬ 
tion [8] shows the strict nronotonicity of the active regions under additional 
assumptions. 


Acknowledgement: We thank Inwon Kim for helpful discussions and in¬ 
terest in this work. 


2. TWO FORMULATIONS OF OPTIMAL PARTIAL TRANSPORT 

In this section we give a more precise statement of the optimal partial 
transport problem that was initiated by Caffarelli and McCann J3}. We will 
provide two equivalent formulations. 


2.1. Notation. Optimal transport problems consist of three basic ingredi¬ 
ents: source and target distributions, cost, and transport plans. We explain 
below the necessary notation for these ingredients we will use in this paper. 
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1. Source and target mass distributions: Let /, g G L 1 (R n ) be two 


nonnegative functions, with 



For simplicity, we will also assume that both / and g have finite second 
moments, namely, f Rn \x\ 2 f(x)dx, f Rn \y\ 2 g{y)dy < oo. 

2. Transport plans: 

a. Let T(f,g) be the set of all Borel measures on R n x R n whose left 
and right marginals are / and g respectively, that is 7 € T(f,g) if 7 
satisfies 



for all A C R n Borel. 

b. Let r <(f,g) be the set of all Borel measures on R n x R n whose left 
and right marginals are dominated by / and g respectively, that is 
7 G r<(/, g) if 7 satisfies 



for all A C R n Borel. 

3. Define now the total mass of a Borel measure 7 on R n x R n by 



4. Costs and modified costs from the transportation cost cap: 

a. Let c : R n x R n -> R be the transportation cost. 

b. FixO < A < 00 . We use this value A as a cap on the value of cost 
that we allow for transportation. This leads us to define the following 
modified transportation cost. 



( 2 . 2 ) 


2.2. Two formulations. Now, we describe the two equivalent formulations 
of the partial transport problem. 

2.2.1. m-problem: changing the portion of mass. The optimal par¬ 
tial transport problem can be stated as follows [3j: and fix 0 < rn < 
minlll/H^i, || 5 || l i} which is the amount of mass to be transported. Given a 
cost function c : R n x R n —>• R the associated cost functional is given by, 
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The partial optimal transport problem consists then to minimize C( 7 ) among 
all possible 7 £ T<(/, g) subject to added constraint of 01 ( 7 ) = m , that is 

(2.3) inf £( 7 ). 

7er < (/,g)& m(7 )=m 

This problem can also be formulated in the following way: Let 
Kf(m) = {h £ L 1 |0 < h < f a.e. & ||h|| L i = m }- 

Then, 


inf C( 7 ) = inf C( 7 ) 

7 gr<(/,g)& 111(7 )=m ')&T(h\,h2)&chi&Kf(m,),h2&K g (rn) 

In particular, when c(x,y) = \\x — y | 2 , the problem on the right corre¬ 
sponds to finding the minimum Wasserstein W 2 distance between the two 
sets Kf(m ) and K g (mn). This point of view is borrowed from a recent work 
of DePhilippis, Meszaros, F. Santambrogio, Velichkov m, where they con¬ 
sidered a projection problem of a probability measure to the set Kf with 
respect to the W 2 metric. 

Existence of the optimizer of (|2.3I) follows easily from a compactness 
argument. Uniqueness of the solution to the optimal partial transport prob¬ 
lem has been established under reasonable conditions. Also, a progress has 
been made in understanding regularity of the free boundaries; , es¬ 
pecially when the cost function c(x, y ) is given by the distance squared 
c(x,y) = 7 dist 2 (x, y). Kitagawa and Pass pTj then considered the problem 
in the case where there are finitely many mass distributions to be partially 
matched, and made a connection to the barycenter problem in the space of 
probability measures that was considered by Agueh and Car her [1]. 


Remark 2.1. We note that the barycenter problem in the space of prob¬ 
ability measures was extended in m to a continuous family of measures 
over Riemannian manifolds, finding Jensen type inequalities over the space 
of probability measures; the corresponding partial transport problem has not 
been considered yet. 


2.2.2. A-problem: changing the transportation cost cap: (pay flat 
rate or do not transfer after a certain threshold). As shown in [3] the 
m-problem (12.31) can be equivalently formulated using a Lagrange multiplier 
A > 0 and considering the unrestricted minimization problem 

(2.4) C\(f,g) = inf f (c(x,y) - X)d-f{x,y). 

Given m, there is A, so that the solutions to these problems coincide [3]. 

We now observe that the value of A carries an economical meaning. First, 
consider 

inf / c x {x,y)d'j(x,y). 

7^r (f,g) 


(2.5) 
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Note that here the admissible set is T(f,g), not T<(/, g). The cost c\ 
reflects the practical situation where for example, a taxi driver charges a flat 
rate after a given upper bound. The following simple observation explains 
the equivalence between (12.41) and (12.51) . thus equivalence between the m- 
problem (12.3|) and the A-problem (12.5p . 


Proposition 2.2 (Equivalence of two A problems). The two problems (I2.4j) 
and m are equivalent, namely, 


( 2 . 6 ) 


7GP 


inf / (c(x, y ) - X)<Ty{x, y) 

c\(x,y)dj'(x,y). 


= —A + inf / 

7'er (f,g) 


Moreover, a minimizer of (12.41) is the restriction of a minimizer of (12.51) on 
the set {(x,y) G R n x M n | c(x,y ) < A}. 


Proof. For each 7 ' G T(/, g), let 7 G T<(f,g) be its restriction to the set 
{(x,y) G I" x R” | c(x,y) < A}. Now, consider 



c\(x,y)d'/(x,y)= c(x, y)dj(x, y) + A / d(n' - 7 )(x, y) 

J I“xl" JR n xR n 

= (c(x,y) - X)d'y(x,y) + A/ d~/(x,y) 

J M"xR" Jl"xl" 

= / (c(x,y) - X)d'y(x,y) + A (from ([24])). 

JR 7 




This shows that the right-hand side of (12.61) is greater than or equal to 
the left-hand side. For the other inequality of (12.61) . notice that for each 
7 G T<(/, g), there is a 7 ' G T(/, g) whose restriction to the set {(x, y) G M n x 
M n | c(x,y) < A} is 7 : in particular, one can take 7 ' = 7 + (/ —/ 7 )<£>(<? —<? 7 ), 
where / 7 , are the marginals of 7 . This completes the proof, including the 
claimed correspondence between the minimizers of (12.41) and (12.51) . □ 


3. Preliminaries, assumptions, and further notation 

In this section, we explain relevant previous results, key assumptions of 
this paper, and further notation. 

The existence and regularity of the optimal plan was done first by Caf- 
farelli and McCann in (3] under the additional hypothesis of the supports 
of / and g to be strictly separated by a plane. In [5],[6] Figalli is able to lift 
the restrictions and even further considers the case when / and g overlap 
on an open set. The strategies used in [3] and [HE] differ greatly. In [3] the 
authors used a Lagrange multiplier A, add a point at infinity and study a 
full transport problem for a redefined cost function. The approach in w 
is to study the convexity problems of the total cost C(m) of the optimizer 
with the given m, and in this way solve the problem directly. Our paper 
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follows the approach of [3], especially because we treat A as a parameter of 
the partial transport problem. 

Assumption 3.1. We now state basic assumptions we impose throughout 
the paper: 

1. (cost function) From now on, let c(x,y) = \\x — y | 2 on W 1 x R n . This 
is the most important special case of the more general class of examples 
(see Remark \3.2\) . We focus on this case for simplicity of treatment. 

2. (domains) fl, A are smooth, open, connected and bounded domains. 

3. (separation of source and target) We assume that there is a hy¬ 
perplane that separates and A. Moreover, we assume that there is a 
constant 6 > 0 that \x — y\ > 5 for all (x, y) € x A. 

4. (measure ~ volume) p = f, u = g, and the functions f,g € L 1 satisfies 
Pi < f,g < P 2 for some constants Pi, @2 > 0 on fl, A, respectively, and 
spt / = 11, spt g = A. 

Remark 3.2. One can generalize most of the results in this paper to the 
cost functions c : R n x R n —> R that satisfy 

(1) (smooth) c € C 2 ; 

(2) (twisted) the maps y i-a \7 x c(x,y) and x i-a V y c(x,y ) are one-to- 
one for all ( x , y) € R n x R"; 

(3) (nondegenerate) the mixed second order derivatives [Df y c(x,y)\ 
give an invertible matrix. 

One exception is the assertion 5 in Theorem S3 where the strict mono¬ 
tonicity of free boundary movement is shown: for this result, we used the 
special structure of the cost function c(x,y) = ^\x — y\ 2 in a essential way, 
to use the theory of classical Monge-Ampere equation. 

The above conditions on the cost function are the usual assumptions one 
require, especially in the regularity theory of optimal transport, though in the 
latter one require more assumptions theory, such as the Ma- Trudinger- Wang 
condition UMJM (see also fWf). All these assumptions are satisfied by the 
quadratic cost c(x,y) = — y\ 2 which is the main focus of many recent 

papers, including this paper. 

3.1. Minimizers, active regions, and potential functions. This sub¬ 
section is devoted to explain some of the fundamental contributions from 
Caffarelli and McCann [3] that are relevant in our paper. As shown in [3], 
under Assumption 13.11 there is a unique minimizer of (12.411 , and for 

(3.1) 7 a := argminCA(/,#), 
its mass satisfies the relation 

/oon ^ 9C\(f, g) 

(3.2) m{ A) := m(7 A ) =-^-, 

moreover, m(A) increases continuously from 0 to min{||/||^i, ||}- There¬ 
fore each mass m can be attained for an appropriate value of A. 
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To understand the behaviour of 7 a, we let, as in m, fx < f,g\ < 9 
denote the marginals of the minimizer r y\ £ of ()2.4ji . Then, the 

regions spt f\ , spt g\ can be regarded as the (closure of the) active regions in 
the source and target domains, where the actual transport occurs. Therefore 
both the source and target domain can be decomposed as an active region 
and an inactive region, and the common boundary arises naturally as a 
free boundary of the partial transport problem. More precisely we have the 
following definition: 

Definition 3.3. 1. (Active regions) The active regions of the partial trans¬ 
port problem are given by the following: 

A \ 1 (x, y) £ spt 7 A }, 

A a : = {v € ^ |3x, (x, y) £ spt 7 a} 

Notice that spt/ a, spt^A are the closures of A^, Arespectively. When 
its meaning is clear from the context, we will simply use Ai to denote the 
active regions. 

2. (Free boundaries of the active regions) Let 

F% = dA% \ dn, F^ = dA t \ dA. 

One of the important observations made in 0 Corollary 2.4] (see also ( 6 j 
Proposition 3.1, Remark 3.3]) is the following interior ball condition for the 
active regions, which under Assumption 13.11 can be stated as 

(3.3) spt f\ = closure of U {x £ LI | c(x,yi ) < c(xi,yi)}, 

(zi,i/i)espt 7 A 

spt g\ = closure of U {y e A I c(xi,y) < c(x i,yi)}. 

(xi,j/i)espt7A 

Another important consequence from (H Corollary 2.4] is that 

(3.4) ./' = fx, 9 = 9\ on spt /a, spt g x , respectively, 

which is important throughout the paper. We also know from [3l Theorem 
3.4] that the active regions increase monotonically along with the parameter 
A and therefore also with rn. This will be an extremely important tool when 
comparing optimal plans associated to different m’s, A’s: 

(3.5) spt / Al C spt fx 2 , spt g Xl C spt gx 2 ; 

for A, < A 2 . 

3.1.1. Augmentation with infinity. In [3] Caffarelli and McCann intro¬ 
duced an effective way to treat the partial transport problem by adding an 
auxiliary point at infinity: Namely, the strategy is to attach a point 60 to 
M n and extend the cost function c 

. I c(x, y) — A if x 7 ^ do and y / c»; 
c(x,y) = < 

[ 0 otherwise, 
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and the measures dy, = f(x)dx and du = g(y)dy to R n = R n U { 6 b} by 


A = T + 11 511L 1 ) 


i > = ||/|| l k5oo 


Then, the problem (12.41) is equivalent to 
(3.6) 



where similarly as before T(y, 0) is the set of Borel measures on R n x R n 
with marginals y, v. In particular, the minimizer 'j\ of (12.41) is nothing but 
the restriction of the minimizer 7 a of (13.61) to R n x R n . We refer [3j Section 
2 ] for more details. 

In particular, from (13.41) this implies that points outside the active re¬ 
gions are matched with 6 b by the plan 7 ^. In other words, the difference 
7 — 7 / is a measure supported on 


({spt(/ - /a) x 6 b}) U ({ 60 } x (spt (5 - g x )) U ({ 60 } x { 6 b}) 
inside the product space R n x R n . This can be written as 
(3.7) 7a - 7a € T< (/ - / A , ||/|| £ i5db) 


+ r< (\\g\\ L ^6b,g - gx) + r< (IbllLi^db, Wfh^db) 


where T<(a,f3) denotes as before the set of measures whose marginals are 
bounded above by the measures a, /3. 

3.1.2. Potential functions. We will make use of the c-potential functions 
in [3]. In our presentation, we will use c or c-convexity instead of c-concavity 
in [3], since convex functions are more natural in terms of Monge-Ampere 
equations that we will consider later. 

Definition 3.4 (c-convexity and c-transform). 

(1) A function u : R n —> (— 00 , 00 ] is said to be c-convex if it is not 
identically 00 on spt (jj) and satisfies 

u(x) = SUp ~c(x,y) - Uc(y) = Ucc, 

y&spt(u) 

where u c is the c-conjugate: 


u c {y) = sup —c(x,y) — u(x). 

x^spt(fi) 

(2) The c-subdifferential of u is defined as 


d c u(x) = {y <E R n | u(x) + u £ (y) = -c(x,y)}, 


d c u = {(x,y) £ I" x 1 " | u(x) + u £ (y) = —c(x,y)}. 

The c-subdifferential d c u c of the c-conjugate u c is similarly defined 
and from symmetry, we have 


d c u = d c u c . 
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Remark 3.5. We obtain definitions for c-concavity as in m Definition 2 . 1 ], 
by just changing the sign and taking supremum instead of infimum above. 

Note that the c-convex functions are locally Lipschitz and semi-convex 
on M n , since they are obtained by taking supremum of C 2 functions (since 
—c(x,y) = —\\x — y \ 2 + A on M n ). In particular, this implies that u are 
differentiable a.e. and they are subdifferentiable, i.e. at each point x £ R n , 
the subdifferential 

du(x) = {p £ 1" | u(y) < u(x) +p ■ (y — x), Vy £ M"} 

exists. This clearly holds also for u c . Recall that, from a well-known theorem 
of Alexandrov, semi-convex functions are twice differentiable a.e. 

The point of the c-convex functions is the following fundamental fact: 

Lemma 3.6 (Uniqueness and characterization). (See J21 Lemma 2.3, Propo¬ 
sition 2.9].) Use Assumption 1 3 . 11 There is a c-convex function u (that is 
unique up to additive constants) such that for each optimizer 7 of (13.61) . 

(3.8) sptf cd £ u{= d £ u £ ). 

Remark 3.7. The above inclusion ()3.8|) can be improved, when restricted 
to SI x A, to 

(3.9) spty fl (fl x A) = d c u n (fl x A), 7 -a.e.. 

The reason is that inside fl x A, d c u is equal to the c-subdifferential d c u of the 
c-potential function u (which is nothing but the restriction of the definition 
of d c u without having the point 60 ), which is well-known (see e.g. w to 
coincide 7 -a.e in fl x A, with the graph of a Borel measurable map T, which 
then coincides with spt 7 - (The map T is the optimal map in the ordinary 
optimal transport problem). 

3.1.3. Normalized potential u\. Since we will consider the optimization 
problem (12.4j) (equivalently (13.61) . for each A, let 7 a, 7 a, denote the minimiz- 
ers of (12.41) . (13.61) . respectively, and u\, the corresponding c-convex function 
for 7 a given in Lemma 13.61 Let u\ denote the restriction of u x to M n , and 
from now on, we use the normalization 

u\(oo) = 0. 

For clarification, we call the restriction u\, the normalized potential. An 
important observation for normalized solutions is the following. Namely, 
from (|3.7I) and (13.81) . we see that 

60 £ d c u c x (y) for each y £ M n \ 

thus, for the normalized solution u\, we have for each y £ M n \ A^, 
u c x (y) = —c(do, y) - u x (ob) = 0 + 0 = 0, 
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using the definition of the cost c and the c-subdifferential. Therefore, we see 
from the definition of d c u\ and u\, 

(3.10) u\(x) > —c(x, y) + A whenever x € M n , and y 6 M n \ 

u\(x) = —c(x,y) + A if moreover y e d c u(x). 

3.1.4. C? or regular partial transport mappings. We finish this section 
with a characterization and regularity result given in [3]: 

Theorem 3.8 (Caffarelli and McCann [3] ; see also mm- Use Assump¬ 
tion E3 Then, for each 0 < A < oo, the optimal transport 7a, i.e. the solu¬ 
tion to m, uniquely exists and is given by a Cjf c mapping T\ : A ^ — > A ^ 
such that 

7a = (id x T a ) # / a 

and 

T\(x) = x + \7u\(x). 

Here, a > 0 depends on n,f,g. Moreover, the free boundaries Fp,Fp are 


4. Main results 

We now explain in detail our main results. First, we introduce some 
notation for simplicity of the presentation. 

Definition 4.1 (Notation <,>,~). We use a > b, a < b, a ~ b, to denote 
a > C\b, a < C 2 &, a > C\b &; a < 6 * 26 , respectively, for some constants 
C\,C 2 , depending only on Q,A and the lower and upper bounds of y, u on 
Tt, A, respectively and in particular, not on the parameters X, m. 

4.1. Main question. The main aim of this work is to get quantitative re¬ 
sults on how the active region changes under the variation of the parameters 
m or A. We give further notation and a few definitions to set up the question 
more precisely, 

Definition 4.2 (Notation and definitions). 1. For i = 1,2, consider the 
values Aj, to* = m(Aj) with 0 < Ai < A 2 0 < m± < m 2 < 1. We let 
7A; 7 i, Ui, v,i, i = 1,2, denote the corresponding objects 7a ; u\,u\ in 

the previous section. By the same way, we let Af, Fp, Fp, denote 
the corresponding active regions and free boundaries. When its meaning 
is clear from the context, we will simply use Ai, Ft. 

2. (b-distance) We define the b-distance between the free boundaries of the 
active regions (in flj: 

Dist b (Fi, F 2 ) := inf{e | Fj C N t (F t U dM), i,j = 1, 2} 

This distance is similar to the Hausdorff distance, but it considers the 
effect from the boundary dQ. 
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Remark 4.3. Note that smallness of Distb (Ti, F 2 ) implies the two free 
boundaries F\ and F 2 are close to each other in a uniform manner outside 
a small neighbourhood of the boundary of the domains. Incorporating the 
boundary dUl for measuring the distance between the free boundaries looks 
somewhat technical, but our method in the proof of Theorem \ f . 5 \ assertion 
3 requires it. However, one can view this distance Distb as the Hausdorjf 
distance between the boundaries of the inactive regions tl\ Ai, i = 1 , 2 . 

We now use the above notation to state the main question of the present 
paper: 

Question 4.4 (Main question). 

• Can we estimate Distb(Fi, F 2 ) in term of the values mi, m 2 or X\, A 2 ? 

• Since clearly the cost cap A and the mass m are related, can we find 
a quantifiable relation between them? 

4.2. Main Results. Addressing Question 14.41 our main results give quan¬ 
titative estimates for the free boundary movement as the parameters m and 
A vary. 

Theorem 4.5 (Quantitative results on free boundaries). Use Assumption \3. 11 
Then, the following hold: 

1. (Lower bound of speed of free boundary movement) 

m 2 -m 1 < Distb(T i,F 2 ); 

2. (Holder continuity of free boundary movement in m) 

Distb (i 7 !, F 2 ) < (m 2 -mi) 7 ; 

3. (Lipschitz bound of free boundary movement in A) 

Distb (Fi, F 2 ) < (Ai — Ai); 

4- (Lipschitz bound of m( A)) 

m 2 - mi < A 2 - Ai. 

5. (Strict monotonicity of free boundary movement) Assume 
further that f,g € C a for some a > 0, then 

Fi(1F 2 = 0. 

Proof. The assertion 1 and 2 will be shown in Section [5j The assertion 3 
will be shown in Section [TJ The assertion 4 follows immediately from the 
assertions 1 and 2. Finally, the assertion 5 will be shown in Section [SJ □ 

An immediate and interesting consequence of Theorem 14.51 Assertion 4 
together with (13.21) is this: 

Corollary 4.6. Use Assumvtion \3.1\ The cost function C\ given by 
C\(f,g)= inf [ [c(x,y) - \]d^(x,y), 

7er <(/,<?) J 

l"xl" 
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is a C 1,1 function as a function of X. 

Remark 4.7. Note as observed in m that the function A 1-4 C\(f,g) is 
concave because it is an infimum of a family of linear functions in X. 

Notice that the Holder exponent 1 jn in the assertion 2 in Theorem 14.51 
is sharp, as seen the example when the target measure is given by (an 
approximation of) a Dirac mass, when m is close to 0. 

Remark 4.8 (Dirac delta target). Consider the case that p, is just the uni¬ 
form density in -Bi(O) and u is Dirac delta concentrated in 2e n . We note 
that the free boundary is always of the form 

F = (u x£A B r{x) (x))nn. 

Therefore, since the target region is just one point we conclude that the 
active region is given by -Bi(O) (~l B r (x). We note from this example several 
remarkable properties. 

(a) The free boundary is strictly monotone. 

(b) The separation of the free boundary can be controlled directly by the 
mass. More precisely we have that the free boundaries are withing 
a tubular neighborhood of size mf! n . Let 0 < m\ < m 2 < 1, then 
the associated free boundaries for the partial optimal transport with 
masses m\ and m 2 are given by 

Dist b (i ? i, F 2 ) < (m 2 - mi) 1/n . 

However, we have the following conjecture: 


Conjecture 4.9. Use the same assumptions and notation as Theorem \f.5\ 
We conjecture the following reinforcement of Theorem f.5. which are all 
related to each other: 


(1) m 2 — m\ ~ A 2 — Ai for | < m; < §; 

(2) Dist b (i ? i, F 2 ) ~ m 2 - mi for ± < m t < § ; 

(3) For Tj < m t < there exists 5 such that 5 > 0 with 6 ~ m 2 — m± 
and Ns(F\) n Ng(F 2 ) = 0. 


This conjecture predicts that the movement of the free boundary and 
the growth of the mass are linear in A. Here, the condition ^ | is 

required as it can be easily seen from the case where v is close to a Dirac 
mass 5 yo and when m is close to 0 or 1: see Remark 14.81 The assertion 5 of 
Theorem 14.51 can be regarded as a partial result for the item (3). 


5. Mass m and the free boundary: Proof of Theorem 14.51 

ASSERTIONS 1 AND 2 

We now start proving the assertions in Theorem 14.51 In this section, we 
focus on the free boundary in the source domain D. Because the assumptions 
we made are symmetrical for D and A, exactly the same proof shows the 
statement for the free boundary in A. 
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We first show the assertions 1 and 2. These follow easily from the prop¬ 
erty (13.31) and simple geometric arguments. First, observe the following fact: 

Lemma 5.1 (Free boundary Lipschitz and semi-convex). Use Assump¬ 
tio ns. 1\ The free boundary F is uniformly Lipschitz and semi-convex. The 
Lipschitz and semi-convexity constants depend only on f l and A. In partic¬ 
ular, / H n ~ 1 (F) is uniformly bounded. 

Proof. It follows immediately from Assumption 13.11 (especially, the items 2 
and 3) and the interior ball property (13..'ll) . The uniform boundedness of 
TL n ~ l {F) is due to Lipschitz property of F and boundedness of II). □ 

Proof of Assertion 1 of Theorem 14.51 From Assumption 13. II item 4, to 
get upper bound on the difference m 2 — mi, we only need to estimate the 
volume difference between the two active regions A\ and A 2 . From the 
Lipschitz property of F\ and F 2 due to Lemma 15.11 the volume difference is 
estimated as 

vol(A 2 \ Ai) < [H n -\Fi)+n n -\F 2 )\ Dist b {F l ,F 2 ) 
which then implies 

m 2 -mi< [n n -\Fi)+'H n - 1 {F 2 )\ Dist b (Fi, F 2 ). 

From the uniform boundedness of TL n ~ l (F) (Lemma 15.ip . we get 

m 2 - mi < Dist b (Fi, F 2 ). 

This proves the assertion 1. □ 

Proof of Assertion 2 of Theorem 14.51 For the two free boundaries F\ 
and F 2 , let x € F 2 be the point realizing the distance d = dist(F\, F 2 ), i.e. 
B c i(x)nFi = 0. Thanks to the interior ball condition (13.31) . then there exists 
y € A j 4 2 such that B dist i x ^(x) n fi C A 2 . This implies 

B d (x) n B dist ( X: y) nfici 2 \4 

Note that 

vol (B d (x) n B dist(Xty) n H) ~ d n , vol(A 2 \ Ai) ~ m 2 - mj. 

Here we used Assumption 13.11 This shows 

Dist b (F 2 , Fi) < (m 2 - mi) 1/n 

as desired. □ 


6. Monotonicity of the potential functions 

For the assertion 3 of Theorem 14.51 we need monotonicity on the po¬ 
tential u\ with respect to A. The goal of this section is to prove such 
monotonicity: 

Theorem 6.1 (Monotonicity). Use Assumption WA 1 Letu\ x andu\ 2 be the 
normalized solutions of the partial optimal transport problem (ESD associated 
to Ai and X 2 , respectively, Then u\ 1 < u\ 2 in Ll. 
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For the sake of a clearer exposition from now on we let 
(6.1) Vi(x) = u Xi (x) + ^|x| 2 . 

An important observation is that vfs are convex for c-convex functions uf s 
(here c(x,y) = \\x — y\ 2 ). Moreover, dvfix) = d c u\ 1 (x), thus we can regard 
Vt>j as the transport map for each partial transport problem (12.51) with 
Ai,A 2 , respectively. Here, each Vu,, i = 1,2, is viewed as a measurable 
mapping defined a.e. on M n . Now, let us recall a version of the Aleksandrov 
Lemma as stated in [Ill- 

Lemma 6.2 (Alexandrov’s lemma). Let v\,V 2 ■ —>• M be convex func¬ 

tions, differentiable at a point p with v\ (p) = V 2 (p), and Viq (p) / Vt >2 (p). 
Define M = {vi > V 2 } and X = {X7vf l (dv\ (M))}. Then X C M, while p 
lies at a positive distance from X. 

We use this lemma to show Theorem 16.11 

Proof of Theorem 16.11 Recall the definition of the convex functions Vi in 
(16.11) and note that it suffices to show v\ < V 2 in fL 

We will proceed by contradiction. Suppose that v\ and V 2 cross in Q, 
i.e. M = {v\ > V 2 } 0 • We hrst assume that they cross in a transversal 

fashion, namely, there is a point p £ such that v\(p) = V 2 (p) and v\,V 2 
are differentiable at p with Vui (p) Vr >2 (p). We note that p € dM, M 
is a subset of the active region of A\ for Ai. Let Y = dv\ (M). Due to 
Lemma 16.21 we know there is an open neighborhood of p, J\f p C D such that 
Xvf l (Y) C M and excludes J\f p . This implies 

M ((Vu 2 )- 1 (T)) < /x(M) < //((Vui)- 1 ^)). 

The hrst inequality comes from the fact that flA/),) > 0 and the neigh¬ 
borhood Af p is not included in Xvf 1 (Y). The second inequality follows from 
the fact that v\ is differentiable a.e., so ji(M \ (Vui)” 1 (y)) = 0. Since 
Y C C Aff we get a contradiction, as both Wi and Vu 2 push forward // 
into v in that region Y. 

In the case that v\ and u 2 cross in D, but does not intersect transversally 
in the interior, we proceed as follows. Lower the function v\ down a bit, i.e. 
consider v\ = v\ — e for each small e > 0 (here, is not normalized anymore) 
and note that the active region and the mass moved by Vwf is the same as for 
v\, because they depend only on Vui. Since the crossing was in the interior, 
then for small e, vf and u 2 will still cross in the interior of fL Consider the 
family {uf} e >o. First, notice that 

u e >oK = M n {v\ > v 2 } = {ui > v 2 }. 

Suppose each set {v\ = u 2 } H {ui > u 2 } has no point where the crossing 
between v\ and u 2 is transversal. Then, it implies that Vrq = Vu 2 a.e. in 
{ui > u 2 }. This means v\ = u 2 in {ui > U 2 }, an obvious contradiction. 




DYNAMICS OF OPTIMAL PARTIAL TRANSPORT 


15 


So, for some e > 0, v\ and V 2 cross transversally in the interior of Cl. Fi¬ 
nally, we can apply the same reasoning as before, to v\ and V 2 , to derive a 
contradiction. □ 


7. Cost cap A and the free boundary: Proof of Theorem 14.51 

ASSERTIONS 3 

In this section, we prove the following proposition, which verifies the 
assertion 3 in Theorem 14.51 

Proposition 7.1. Use Assumvtion 17.11 Recall the free boundaries to 
the solutions of (12.41) . with Ai < A 2; respectively. Then 

(7-1) F 2 A C AT c(A 2 _ Al) (FiUdA), 

Ti A ClV c(A 2 _ Al) (F 2 UiM), 

for some constant C, depending only on Cl, A. Note that this implies 

Dist b (F A ,F A )<A 2 -A 1 . 

Since Assumvtion \3. 1\ is symmetrical, we also have 

Dist b (F 1 n ,if)<A 2 -Ai. 

Proof. Let u\ and U 2 be the normalized potential functions to the solutions 
of (12.41) . with Ai < A 2 respectively. 

We begin by showing the first inclusion in ( 17 . 11 ) . Consider an arbitrary 
point 2/2 G F <2 ■ It has a corresponding point X 2 G clST (in fact in dCl due 
to [3, Corollary 6.9]), such that y 2 G d c U 2 (x 2 ). Consider the line X 2 IJ 2 , and 
then either x 2 ?/2 H T A = 0 or there exists y\ G X2V2 D F\ A . In the first 
case, there is a point in X 22/2 H dA D (M n \ A a ), and we denote it also by 
y\. Notice that the functions m,i(x ) = ~\\x — yi \ 2 + A i, i = 1,2, satisfy 
mi < Ui, and m2(2:2) = 112(^2), since we can apply (| 3 . 10 j) to 2/2 G d c ii2(x2) 
and yi GC M n \ A a . Therefore, noting that u\ < 112 in Cl from Theorem 16.11 
—c(x 2 ,y 2 ) + A 2 > —c(x 2 ,yi) + Ai, which gives us 

c(x 2 ,y2 ) - c(x 2 ,yi) < A 2 - Ai. 

Since c(x,y) = \\x — y \ 2 and II and A are separated (3 of Assumption 13.II) . 
using the fact the points X 2 ,yi,y 2 are on the line T 2 I/ 2 , we have then 

c(x 2 ,y 2 ) - c(x 2 ,yi) > I2/2 - yi|, 

proving the desired result F ^ C A r c -( A 2 _ Al )(T 1 A U dCl). 

By interchanging the role of F and T A , the exactly same method shows 
the other inclusion F A C N c ^ 2 _ Xl ' ) {F^ U dCl), completing the proof. □ 
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8. Strict monotonicity of free boundaries: Proof of 
Theorem 14. 5 L assertion 5 


A desired result on the free boundary movement is to get a quantitative 
separation of free boundaries, showing points in two different free boundaries 
corresponding two different values of rri or A, are away from each other 
by a positive distance whose lower bound is controlled by m or A. As a 
partial progress toward this direction, we show in this last section, the free 
boundaries indeed do separate from each other, but without a quantitative 
estimate, which also establishes the strictness of the monotonicity of the 
active regions. To do this we require the following technical assumptions, 
which are mainly for using the established regularity theory of the optimal 
partial transport [3]: 

Assumption 8.1 (Additional assumptions for strict monotonicity). 

1. (smooth and convex domains) 12, A are strictly convex and bounded 
domains with smooth boundaries. 

2. (smooth densities) The functions f,g £ C 1 . 

For /, g £ C“ c , consider v(x) = u(x) + \\x \ 2 , where u is the potential 
function to the solution of the partial transport problem (12.411 . Then, v 
satisfies in the interior of the active region, det(Z? 2 u) = f(x)/g(Vv ),and v 
(thus u ) is C 2 0 ’“; this interior C 2 regularity is due to Caffarelli [2|. The reason 
we require a further regularity C 1 is to apply the Hopf lemma (Lemma 18.2D 
and the strong maximum principle (Lemma 18.3D shown below using the C 1 
assumption on / and g. 


8.1. Hopf’s lemma and the strong maximum principle. In this sub¬ 
section for reader’s convenience, we give a short proof of a Hopf Lemma and 
strong Maximum Principle for the Monge-Ampere equation 

(8.1) det (D 2 v(x)) = F(x,Vv(x)), 

which should be known to experts. We assume that F is C 1 , which for our 
partial transport problem correspond to the case where the densities f,g of 
the measures /r, u, respectively, are C l , and g is bounded from below. 


Lemma 8.2 (Hopf lemma for Monge-Ampere). Let B be an open ball and 
vi,V 2 £ C 2 (B) n C l (B) be uniformly convex solutions of in B, where 
F is a differentiable function. Suppose V 2 > v\ on B and there is a point 
xq £ dB such that ui(xo) = 112 ( 2 : 0 ) • Then, we have the strict inequality 


d y 2, w dvi . ■ 

where ft is the inner normal vector of dB and xq . 


Lemma 8.3 (Strong maximum principle for Monge-Ampere). Let U be an 
open connected domain and V\,V 2 £ C 2 (U) be uniformly convex solutions of 
(EH in U, where F is a differentiable function. Suppose V 2 > v\ on U and 
there is a point xq £ U such that vi(xq) = V 2 {xq), then v\ = V 2 on U. 
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Proof of Lemmata 18.21 and 18. 3L Following a standard argument below, 
we reduce to the classical Hopf’s lemma and the strong maximum princi¬ 
ple of the liner uniformly elliptic equations. First we note that since the 
solutions are C 2 , then the Monge-Ampere equations are uniformly convex 
(non-degenerate). We observe that for w = vi — v 2 , 

0 = det(Z) 2 ui) — det(-D 2 u 2 ) — (F(x, Vrq(x)) — F(x, Vv 2 (x))) 

f 1 d 

= J — [det((L> 2 (fui + (1 -t)v 2 )(x)) + F(x,V(( 1 - t)v i + iu 2 ))(x))] dt 

= a,ijD 2 jW + CjDjW 

where 

a%j = [det D 2 (tvi + (1 - t)v 2 )\ [D 2 (tv 1 + (1 - t)v 2 )]^ 

Ci(x) = [D Pi F] (x,V(l -t)vi(x) +tv 2 (x)). 

Here, D Pi F is the derivative of F(x,p) in the p variable. We notice that the 
coefficients a^j are uniformly elliptic U due to uniform convexity and the C 2 
regularity of v\, v 2 . Now, one can use the classical Hopf’s lemma and the 
strong maximum principle for liner elliptic equations (see e.g. [ 8 ], Theorem 
3.5]). □ 

8 . 2 . Strict monotonicity of the active regions. Suppose we have two 
Lagrange multipliers X 2 > X±. We would like to study the strict monotonic¬ 
ity of the active regions. We already know that they are all monotone, but 
at this time we can’t discard points for which the free boundary for different 
Lagrange multipliers remain fixed. 

To achieve this we approach the problem in two steps. The first step 
is to note that if we have a point that belongs to both free boundaries, 
then the image through each optimal map is the same; here, we use the 
C l regularity of u along the free boundary from [3] . The second step is 
to use the monotonicity of the solutions plus a Hopf Lemma argument to 
contradict this fact. For the discussion of this section, recall the notation 
Vi(x) = Ui(x ) + \\x\ 2 . 

Lemma 8.4 (Touching implies the same gradient). Under Assumptions\3fJ\ 
and \8.1i let X 2 > X\ and suppose that there is xq such that xq £ dF\ and 
xq £ dF 2 . Then we have Vvi(xq) = Vv 2 (xq). 

Proof. Note that from Theorem 13.81 we have that the optimal partial trans¬ 
port maps are given by F t = Vt>j. Now, let us first consider the case 

V^i(^o) , FIv 2 {x q) 

|Vui(x 0 )| ^ |Vu 2 (* 0 )r 

Then we have two supporting balls with different centers and radii, namely 
B ri (Vv i(xo)) and B r2 (Vv 2 (xo)), whose boundaries cross at xq transversally. 
Now, note that a graph satisfying the previous statement cannot be C l , 
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since no C 1 graph can have two supporting balls. Therefore, since the free 
boundaries are C* 1 , as shown in |3j, this leads to a contradiction. 

Now suppose we are in the second scenario 

Vti(x 0 ) _ Vt 2 (x 0 ) 

|Vi;iOo)| |Vn 2 (x 0 )| 

Note that this implies that the image Vui(xo) and Vu 2 (xo) ni the same 
line from the point xq. But, notice that from [3j Corollary 6.9], they are both 
on the boundary <9A. Thus, strict convexity of A implies Vui(xo) = Vu 2 (xo), 
completing the proof. □ 

We now prove the desired strict monotonicity of the active regions, by 
showing no intersection of the free boundaries occurs for different values of 
A. 

Theorem 8.5 (Strict monotonicity of the active regions). Under Assump- 
tions \3. 1 1 and \ 8. 11 let A 2 > Ai. Then we have that F 2 n F\ fl R = 0. 

Proof. First, we suppose that there is an interior point in the free boundary 
that belongs to both free boundaries, that is, there is xo G F 2 fl F\ n Q. 
Due to Lemma EH we know that Vu 2 (xo) = Vtr(xo). Now, due to the 
monotonicity of u 2 and u\ lTheorem l 6 .ll) we can apply the strong maximum 
principle (Lemma 18.31) in any region in A\ C i 2 , where both functions are 
C 2 , showing that u 2 > u\ (thus u 2 > iq) in A\. Finally, notice that at the 
point xo, the free boundary F\ satisfy the interior ball condition (13.31) . thus 
we can apply Hopf Lemma (Lemma 18.2p to derive a contradiction. This 
completes the proof. □ 
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